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Abstract 

Using observations from the Multi-angle Imaging Spectroradiometer (MISR), the Moderate 
Resolution Imaging Spectroradiometer (MODIS), and the Clouds and the Earth's Radiant Energy 
System (CERES) instruments onboard the Terra satellite; we present a new technique for studying 
longwave (LW) radiative forcing of dust aerosols over the Saharan desert for cloud-free conditions. The 
monthly-mean LW forcing for September 2000 is 7 Wm' 2 and the LW forcing efficiency 1 (LW e ff) is 15 
Wm' 2 . Using radiative transfer calculations, we also show that the vertical distribution of aerosols and 
water vapor are critical to the understanding of dust aerosol forcing. Using well calibrated, spatially and 
temporally collocated data sets, we have combined the strengths of three sensors from the same satellite 
to quantify the LW radiative forcing, and show that dust aerosols have a "warming" effect over the 
Saharan desert that will counteract the shortwave “cooling effect” of aerosols. 


LW tff = Longwave forcing per unit aerosol optical thickness at 0.55 pm 
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1 . Introduction 

Desert dust is considered to be one of the major sources of tropospheric aerosol loading, and play 
an important role in climate forcing studies [Kaufman et al., 2002; Christopher and Zhang, 2002]. 
Widely prevalent in the tropics [Prospero, 1999], dust aerosols are effective in reflecting solar energy 
back to space thereby "cooling" the earth’s surface [Tegen et al., 1996; Myhre et al., 2003]. Besides their 
radiative impact in the shortwave portion of the electro-magnetic spectrum, dust aerosols also have an 
important radiative effect in the longwave (LW) [Sokolik et al., 1998; Liao and Seinfeld, 1998]. Having 
mean particle sizes on the order of several micrometers [Kaufman et al., 2001], dust aerosols can 
effectively reduce the earth’s LW emission by re-emitting at a colder temperature when compared to the 
surface and thereby "warming" the earth [Sokolik et al., 1998; Tegen et al., 1996], However, numerical 
simulations of dust LW forcing remain a challenging task because of the high spatio-temporal variation 
of dust properties [Myhre et al., 2003]. 

Using satellite observations, several studies have attempted to estimate the LW radiative forcing 
of dust aerosols. For example, using broadband measurements from the Earth Radiation Budget 
Experiment (ERBE), the changes in LW and shortwave (SW) radiation at the top of atmosphere (TO A) 
with and without dust aerosols has been reported [e.g. Ackerman and Chung, 1992; Hsu et al., 2000]. 
Since the ERBE instruments were not designed to detect aerosols, narrowband satellite observations 
from the Advanced Very High Resolution Radiometer (AVHRR) [Ackerman and Chung, 1992] and the 
Total Ozone Mapping Spectrometer (TOMS) [Hsu et al., 2000] were used to identify the presence of 
dust plumes within the large ERBE scanner footprints. Although the AVHRR and the ERBE scanner 
were on the same satellite (e.g. NOAA-9), accurate detection of dust aerosols over bright targets such as 
the Saharan desert was not possible due to the limitations of the AVHRR. The TOMS aerosol index 
which is the difference in aerosol absorptivity between two UV channels has been used to detect dust 
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aerosols over bright targets such as deserts, but the technique is insensitive to dust layers at low altitudes 
[Herman et al., 1997] and since the TOMS and ERBE instruments were not on the same satellite, 
precise temporal collocation was not possible. 

With the launch of NASA’s Terra satellite, well-calibrated measurements from the same satellite 
are now available for examining the role of aerosols on climate. For example, using MODIS and CERES 
data, aerosol radiative forcing over the global oceans has been completed [Christopher and Zhang, 

2002]. One of the primary roles of the Moderate Resolution Imaging Spectroradiometer (MODIS) and 
the Multi-Angle Imaging Spectroradiometer (MISR) is to detect aerosols and retrieve their properties. 
The MODIS with 36 channels has improved spectral, spatial and radiometric resolutions when compared 
with previous imagers [King et al., 1992] and has been used to retrieve aerosol properties over the global 
oceans [Remer et al., 2002] and over land areas with low surface reflectances [Chu et al., 2002]. The 
MISR with multi-angle measurements and four spectral channels has also been used to identify aerosols 
and retrieve aerosol properties over the global oceans and land including bright targets such as deserts 
[Kahn et al., 2001]. In this study, we combine the strengths of three instruments on the Terra satellite to 
examine cloud-free longwave dust aerosol forcing. We use the MODIS to identify clouds because of its 
multi-spectral capabilities and higher spatial resolution when compared to MISR and CERES. Although 
the MODIS provides aerosol properties over the global oceans, aerosol retrievals over highly-reflective 
surfaces such as deserts are not available. Therefore we use the multi-angle capabilities of the MISR 
instrument to obtain aerosol optical depth. Finally we use the CERES to obtain broadband TOA 
longwave fluxes. 
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2. Data and Methods 

Four data sets are used in this study; the MISR Level 2 (MIL2ASAE) 17.6 km- daily aerosol product 
[Diner et al., 2001]; the 5 km- MODIS Level 2 daily (MOD06) cloud product [Ackerman et al., 1998] 
and the 30 km (nadir view) pixel level CERES data [Wielicki et al. 1996]. The MISR aerosol product 
contains aerosol optical thickness and other related properties at 0.558pm. The expected uncertainties in 
MISR aerosol retrievals are ±0.05 for aerosol optical thickness (AOT) less than 0.5, and ±10% for AOT 
greater that 0.5 [ Martonchik et al., 1998], Recently, Diner et al., [2001] showed that there was excellent 
agreement between MISR retrievals and ground based sun-photometer measurements over southern 
Africa during August and September 2000. The CERES scanner on Terra observes SW and LW 
radiances that are converted to fluxes using angular distribution models, which take into account the bi- 
directional characteristics of a reflecting surface [ Wielicki et al., 1996]. Spatial collocation is performed 
using point spread functions of the CERES scanner [ Christopher and Zhang, 2002]. 

The area of study is the Saharan desert (15°- 40°N, 20°W- 40°E). The CERES pixels that are labeled 
as “clear desert” by the CERES data were first selected. However, due to the large footprint of the 
CERES scanner, these pixels could still have some cloud contamination [Christopher and Zhang, 2002]. 
To eliminate these cloud effects, the MOD06 cloud product is collocated with the CERES data and only 
those CERES pixels that are at least 99.9% cloud free as identified by the MOD06 data were used. For 
these cloud-free CERES pixels, the best-fit aerosol optical thickness (T0.55) from the MISR data is 
obtained. As noted in Christopher and Zhang [2002], stringent cloud rejection criteria such as this could 
eliminate some aerosol plumes that have high optical thickness. 
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3. Results 

3.1 Satellite Observations 

Figure la shows the monthly-mean MISR T0.55 for cloud-free CERES pixels. Data is gridded every 
l°xl° and smoothed for display purposes. The major dust plumes are clustered in three regions; within 
Mali and Mauritania, Chad and Niger, and Sudan. Areas shown in white are due to missing data, data 
gaps and also due to cloud-screening criteria 2 . Superimposed on Figure la are six selected regions with 
optically thick aerosol plumes. The latitude and longitude boundaries of the six regions are shown in 
Table 1 . For reference purposes, we also compared the spatial distribution of dust aerosols from MISR 
with those from the TOMS [Herman et al., 1997]. Figure lb shows the monthly-mean Earth Probe 
TOMS aerosol index (version 8, level 2) for the same study period. Pixels with reflectivity greater than 
20% in the TOMS 360 nm channel are excluded to avoid possible cloud contamination. The spatial 
distribution of dust aerosols derived from these two independent instruments are consistent over most of 
the study regions. However, over region 4, thick dust plumes shown in figure la are not seen in figure 
lb. Since the TOMS aerosol index is insensitive to dust plumes at low altitudes [Herman et al., 1997], it 
is possible that these aerosols are not detected by TOMS although we have no ancillary information to 
verify this. 

Figure lc shows the MISR optical thickness versus the CERES LW flux for the six selected regions. 
The CERES fluxes are binned for every 0. 1 MISR optical thickness intervals. The thick black line in 
figure 1 shows the mean value for all regions. With the exception of region 4, as MISR AOT increases 
the CERES LW flux decreases because dust aerosols emit at a colder temperature when compared to the 
surface. The reason for the weak relationship between MISR xo.55 and CERES LW flux in region 4 is 
due to the combined effect of dust aerosol vertical distribution, column water vapor, and surface 

2 Data from September 8,17, 24, 27, and 30 was not used due to these criteria. 
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characteristics (see section 3.2). We use this relationship between MISR T0.55 and CERES LW flux and 
calculate dust aerosol forcing efficiency for the six selected regions (Table 1). The dust aerosol forcing 
efficiency (LW e ff) is defined as the LW flux perturbation due to dust aerosols per unit optical thickness, 
and is derived from the difference in CERES LW flux for MISR To. 55 values between 0 and 1. The 
largest LW e ff of 20 Wm' 2 is found over region 2, and the lowest LW e ff of near 0 Wm' 2 is found over 
region 4. The other four regions have values ranging from 11 to 19 Wm' 2 . The variability of LW e ff is due 
to differences in surface and atmospheric conditions (see section 3.2). Averaged over all pixels, the 
regional mean dust aerosol forcing efficiency is 15 Wm' 2 . 

To examine the spatial distribution of the aerosol radiative effect, we define the TOA dust LW 
forcing as F c]r - F a [ Christopher and Zhang, 2002], where F c i r is the TOA LW flux observed in cloud and 
aerosol free conditions, and F a is the TOA LW flux observed over dust regions. The F c i r is obtained by 
averaging the CERES cloud and aerosol free pixels with MISR T0.55 <0.1. Due to the persistence of dust 
aerosols, only one F c i r value is used for each one of the six selected regions. However for the rest of the 
study area, F c i r is obtained by averaging the LW fluxes for each 3°x3° grid and F a values are reported for 
every l°xl° grid. Figure Id shows the dust aerosol LW forcing over the Saharan desert. In general, as 
AOT increases the LW forcing also increases. However, over region 4, 5, and 6, although dust aerosol 
loading is high, the aerosol induced LW perturbations are low. The averaged dust aerosol LW direct 
forcing for the six regions range from 0 to 14 Wm' 2 (Table 1), and the averaged dust aerosol LW direct 
forcing over the entire study area is 7 Wm' 2 . It is interesting to note that over region 4, where MISR 
observed optically thick aerosols that was not observed by TOMS, the averaged dust LW direct forcing 
is near 0. It is possible that the dust plume in this region is at a very low altitude and emits at a similar 
temperature as the surface. 
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3.2 Sensitivity studies 

Although our analysis clearly shows a longwave “warming” effect, the differences in LW flux as 
a function of To.ss among the six regions are due to both atmospheric and surface properties including 
surface temperature (T s ), and vertical distribution of water vapor and aerosols [Hsu et al., 2000; 
Ackerman and Chung, 1992]. In this section, we tested the sensitivity of LW e fr to changes in T s , column 
water vapor, and aerosol vertical distribution over cloud-free desert conditions. A four-stream radiative 
transfer model [Fu and Liou, 1993] was used to compute TO A LW fluxes for various atmospheric and 
surface conditions with desert aerosol properties as described by Hess et al., [1998]. The vertical 
distribution of dust aerosol extinction is assumed to be exponential, and the dust extinction scale height 
is used to represent the vertical distribution. For the atmospheric temperature and water vapor profiles, 
we use rawinsonde measurements [27.20°N, 2.47°E] from September 1 at 12 UTC (hereafter referred to 
as profile 1). The T s is 3 16K and the column water vapor is 1 .6 cm. To test the sensitivity of T s , we used 
another case that has the same water vapor distribution as profile 1, but with different temperature 
profiles from the surface to 700 mb, The T s of the new profile (hereafter referred as profile 2) is 300K. 
We further varied the water vapor amount from 1.6 to 3.2 cm and aerosol extinction scale height from 1 
to 2 km. 

Table 2 lists the LW e ff as functions of aerosol extinction scale height, T s and column water vapor. 
The profiles 1 and 2 are identified by their surface temperatures of 316K and 300K respectively. When 
the dust height is changed from 1 to 2 km while keeping the surface temperature (316K) and the water 
vapor loading the same (1.6 cm) , the LW e ff increases from 1 1 .4 to 18.7 Wm' 2 , because the dust layer 
now emits at a colder temperature and consequently, water vapor absorption above the dust layer is 
reduced. Increasing the water vapor loading from 1.6 to 3.2 cm for the same temperature (316K) and 
aerosol extinction height (1km) reduces the LW e ff from 1 1.4 to 8.0 Wm' 2 because the absorption due to 
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water vapor above the dust layer has increased. Increasing the surface temperature from 300 to 316 K 
for the same water vapor loading ( 1 . 6 cm) and aerosol scale height ( 1 km) will increase the LW e ff from 
6.3 to 1 1.4 Wm ' 2 because the cloud and aerosol free emission temperature increases. We also note that 
the TOA longwave flux is very sensitive to column water vapor. A one centimeter increase in column 
water vapor has a similar effect as a unit aerosol optical depth increase in reducing the TOA LW flux 
[Hsu et al., 2000]. Therefore, the uncertainties in column water vapor could induce uncertainties in the 
derived F c i r . In summary, LW e n is proportional to T s and aerosol extinction scale height, and inversely 
proportional to column water vapor. All of these parameters are critical to the understanding of LW dust 
direct forcing, and precisely collocated observations of the three parameters from other sources are 
needed for future dust forcing studies that utilize broadband satellite observations. 

4. Summary 

In this paper, we have demonstrated a new technique for estimating the dust aerosol LW direct 
forcing over the Saharan desert. Through the synergistic use of MODIS, MISR and CERES instruments 
onboard the same satellite, our study clearly shows that there is a strong dust warming effect of 7 Wm ’ 2 
over the cloud free Sahara desert regions for September 2000 that can offset the shortwave cooling 
effect of aerosols. The regional mean dust LW direct forcing efficiency is 15 Wm ' 2 per T 0 . 55 . However, 
since water vapor dominates the radiative processes in the longwave part of the electromagnetic 
spectrum, simultaneous measurements of both aerosol height and water vapor are needed to reduce the 
uncertainties in the longwave radiative forcing calculations of dust aerosols. 
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Figure Captions 

Figure la. Spatial distribution of MISR optical thickness at 0.55pm over the area of study for 
September 2000. Missing data is shown in white and six selected regions (see Table 1) are shown in red 
boxes. 

Figure lb. Spatial distribution of the TOMS aerosol index for September 2000. 

Figure lc. MISR to. 55 nm vs. CERES LW fluxes for six selected regions. Solid black line shows the mean 
values for the whole region. 

Figure Id. Spatial distribution of CERES -derived LW direct dust forcing over cloud free desert regions. 
Positive values denote a “warming” effect. 
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Table 1. Summary of MISR optical thickness and CERES fluxes for the six selected regions. 


Region 

Latitude 

degrees 

Longitude 

Degrees 

Dust LW 
forcing 

Wm' 2 

Dust LW direct 
forcing efficiency 
Wm" 2 per AOT 

AOT 

* 

(0.55pm) 

1 

16N-28N 

16W-4E 

14 

19 

0.69 

2 

15N - 20N 

5E-22E 

12 

20 

0.75 

3 

15N - 22N 

22E -36E 

9 

19 

0.66 

4 

23N - 32N 

25E - 35E 

0 

0 

0.52 

5 

27N - 33N 

15E-25E 

4 

11 

0.52 

6 

23N - 27N 

15E-25E 

3 

14 

0.44 

Whole area 

15N -40N 

20W-40E 

7 

15 

0.55 


* Different methods are used in deriving dust LW forcing and LW forcing efficiency (see text). 
Therefore, the averaged dust LW forcing divided by the averaged AOT may not exactly match the dust 
LW direct forcing efficiency. 


Table 2. Dust LW direct forcing efficiency (LW e ff)* as functions of surface temperature, aerosol 


extinction scale height and column water vapor content from radiative transfer calculations. 

Surface 

Aerosol extinction 

Column 

LW eff 

temperature 

scale height 

Water vapor 


K 

km. 

cm. 

Wm' 2 per AOT 

316 

1 

1.6 

11.4 

316 

1 

3.2 

8.0 

316 

2 

1.6 

18.7 

316 

2 

3.2 

14.2 

300 

1 

1.6 

6.3 

300 

1 

3.2 

4.2 

300 

2 

1.6 

11.6 

300 

2 

3.2 

9.0 


* Defined as the difference in TOA LW fluxes for AOT values between 0 and 1 . 
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